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1. Introduction 

In order to study the chiral limit in lattice QCD with Wilson Fermions, it is essential to un- 
derstand the phase structure at nonzero lattice spacing. Since the continuum limit, a — > 0, and the 
chiral limit, m — > 0, do not commute, lattice spacing effects can induce new phase transitions which 
have no analogue in the continuum. These new phase structures are known as the Aoki phase [[[]] 
and the Sharpe-Singleton scenario [[2p. The presence of new phase structures due to the lattice arti- 
facts may at first appear to be a severe drawback for the application of Wilson fermions, however, 
these lattice artifacts can be turned to one's favor: The Aoki phase is reached through a second 
order phase transition and just outside the Aoki phase the pions have dispersion relations which 
resemble those of almost massless quarks in the continuum. 

Because the Wilson term breaks chiral symmetry the effects of the lattice spacing lead to new 
terms in chiral perturbation theory. This extended low energy theory is know as Wilson chiral 
perturbation theory. The precise form of the terms was worked out in [Q, f|, [|] and at order a 2 
involves three new low energy constants, which encode the severity of chiral symmetry breaking 
by the lattice artifacts. Recently a large number of new analytic results concerning chiral dynamics 
with Wilson fermions have been obtained from Wilson chiral perturbation theory by working from 
the perspective of the Wilson Dirac eigenvalues [§, [7j g |, ^ [II], [12|, 0, |14[ [B], [l|, [n|]. The 



aim in this review is to give an introduction to the new understanding which the results has brought 
about rather than on the technical details. 

We will focus on three of the most important new insights which follow from the results 
namely, 1) that the Sharpe-Singleton scenario is only realized in unquenched simulations, 2) that 
the a I y/V scaling of the smallest eigenvalues of the Hermitian Wilson Dirac operator, observed on 
the lattice [18,1^, ^], is a good sign, and finally 3) that the results offer a new way to measure the 



additional low energy constants of Wilson chiral perturbation theory. 

Two new results are also included in this review. First, we determine the meta-stable region 
around the Sharpe-Singleton 1st order phase transition, and second we show that the sampling of 
the individual sectors is essential in order to determine the vacuum structure in the Aoki phase. 



2. Wilson Chiral Perturbation Theory in the £ -regime 

The Wilson term in the Wilson Dirac operator 

CLV 

£V = -7m(V m + V*)- t V m V; (2.1) 

breaks chiral symmetry. Therefore, the shortest length scale a also affects the physics of longest 
length scale \/m n . This effect is captured by Wilson chiral perturbation theory [^, |2|, ^ §], for 
which we now give a focused review. 

As always in chiral perturbation theory it is essential to set up a counting scheme. Here we 
will work in the £ -regime where the V — > °° limit is tied to the chiral and continuum limit such that 

mV, £V, and a 2 V (2.2) 
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are kept fixed and of order unity (£ is the axial mass). The partition function of Wilson chiral 



perturbation theory at leading order then reduces to a group integral ||21 



Z Nf (m,C;a)= / dU e s[u \ (2.3) 

JsU(Nf) 



where the action is 



S = ^LVTr(U + + |zVTr(t/ - U^) (2.4) 

-a 2 VW 6 [Tr (U + U^)] 2 - a 2 VW 7 [Tr (U - t/ 1 ")] 2 - a 2 V W 8 Tr(£/ 2 + U^ 2 ). 

The convention for the low energy constants We, Wj and Wg is that of ^ p|. In [|5|] these constants 
are denoted respectively by —We, —Wy' and — W% . 

As in the continuum [23] it is most useful to consider Wilson chiral perturbation theory in 
sectors with fixed index v, such that [0, |8|] 



ZUm,£;a)= [ dU det v V e s[u] , (2.5) 

; JU{Nf) 

where we obviously have 

Z Nf (m,£;a) = ^Z^(m,£;a). (2.6) 

V 

The index defined in the above decomposition of the partition function is the index of the 
Wilson Dirac operator [||, ||] 

v = £sign«*|y 5 |*}), (2.7) 

k 

where \k) is the eigenvector associated with the &'th real mode of the Wilson Dirac operator in the 
physical branch. We will make explicit use of this index in section [7] where we discuss the micro- 
scopic spectral density. In the first sections we will work at mean field level where the dependence 
on the index is suppressed. 

From the action in Eq. (2.4) it is clear that the partition function in the e-regime of Wilson 
chiral perturbation theory only depends on the dimensionless scaling variables 

m = mVL, C = £V£ and a 2 = a 2 VWj. (2.8) 

The mean field limit corresponds to the saddle point approximation to the group integral in the limit 
where these dimensionless numbers are all much larger than unity. 



3. Aoki phase and Sharpe- Singleton scenario 

The form of the order a 2 terms in Wilson chiral perturbation theory are uniquely determined 
by the way in which the Wilson term breaks chiral symmetry. Each of the three new terms comes 
with a new low energy constant, the value of which is not fixed by the flavor symmetries. The sign 
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Figure 1: The mean field potential for the orientation of the vacuum. When cos(0) — ±1 the lattice 
simulation is in the normal phase with the standard chiral condensate. Left W% + 2Wg > 0: The Aoki phase 
develops when the minimum of the potential is between -1 and 1. Right W% + 2W(, < 0: For the Sharpe- 
Singleton scenario the minimum is always at -1 or 1. The additional meta-stable minimum at cos(0) = 
— 1 (red curve) becomes degenerate with the minimum at cos(0) = 1 for m = (blue curve) resulting in 
the 1st order Sharpe-Singleton phase transition. The parameters chosen for the plots are m = 10,5,0 and 
S(al + 2al) = ±4. 



of these three additional low energy constants, however, is determined by the 75-Hermiticity of the 
Wilson Dirac operator 

y 5 D w y 5 =D ] w . (3.1) 

Only for 

W 6 < 0, W 7 < and W 8 > (3.2) 

does the Wilson chiral Lagrangian describe lattice QCD with a 75-Hermitian Wislon Dirac operator 
[Q, ||, 24, 14]. As we shall now demonstrate these signs give vital information on how the Aoki 



phase and the Sharpe-Singleton scenario are realized. To show this let us first work out the phase 
structure of lattice QCD with Wilson fermions at small m and a using Wilson chiral perturbation 
theory at mean field level. With the mean field ansatz for the orientation of the Goldstone field [||] 
(here and below we focus on the case of two mass degenerate flavors) 

U = cos(e) + /sin(0)a 3 (3.3) 

we find the mean field action 

S M f = 2mcos(0) -4^(2cos 2 (0) - 1) - 16^cos 2 (0). (3.4) 

This effective potential for cos(0) is plotted in figure |l| for the two cases W% + 2W(, > (l.h.s.) and 
W$ + 2Wd < (r.h.s.). In the first case the orientation, as given by cos(0), continuously moves 
from 1 to -1 as a function of the quark mass. This is the Aoki phase In the second case, where 
W% + 2W(, < 0, the orientation jumps from 1 to -1 at m = 0. This is the first order Sharpe-Singleton 
phase transition ph. 
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m 




Figure 2: The phase diagram for lattice QCD with Wilson fermions in the (8(Wg + 2W6)a 2 ,m)-plane. 
For positive values of %cr(W% +2W(,) the Aoki phase dominates with the 2nd order phase transition at 
\m\ = 8a 2 (Wg +2We) and associated Goldstone modes for smaller values of \m\. With negative values of 
8a 2 (W8 +2W(,) the Sharpe-Singleton scenario dominates with the 1st order phase transition at m = and 
spinodal lines at \m\ = |8a 2 (W8 +2Wg)|. The meta-stable region is hence between the dashed blue lines. 



The pion masses for \m\L > &(Wg + 2W(,)a 2 are given by the order-a 2 corrected GOR-relation 



,25] 



m 2 F 2 

-^ = \m\L-8{Wg + 2W (> )a 2 . (3.5) 

Approaching the Aoki-phase (where Wg + 2W$ > 0) the squared mass of the pions goes to zero 
linearly with the quark mass. If we absorb the order a 2 -correction into the quark mass this behavior 
of the pion masses is just like in the continuum. On the contrary in the Sharpe-Singleton scenario 
(where Wg + 2W(, < 0) the pion mass is always positive. See the rightmost panels of figure || and 
figure |5| respectively. 

Since the Sharpe-Singleton scenario involves a first order phase transition it is natural to ask: 



What is the meta-stable region? The answer is readily seen from Eq. ( p.4| ) or graphically from 
the r.h.s. of figure [[]. The spinodal line occurs when the maximum of the effective potential enters 
between —1 and 1. That is, at exactly the same quark mass as the Aoki phase would have occurred 
for the opposite sign of Wg + 2W(>. This is illustrated in figure ^[ 



4. A puzzle 

The above analysis shows that either the Aoki phase or the Sharpe-Singleton scenario will 
dominate if we take the chiral limit prior to the continuum limit. Which of the two is realized, 
depends on whether Wg is more positive than 2W(, is negative and hence on the details of the specific 
lattice simulation in question. This analysis was extended to the quenched case in [ 26 ] where it was 
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Figure 3: The physical branch of the spectral density of the Wilson Dirac operator (left) the chiral conden- 
sate (middle) and the pion mass (right). The vertical dashed lines mark the support of the eigenvalue density 
of the Wilson Dirac operator. The location is at |/?e[z]|£ = %(W% + 2W(,)a 2 . The Aoki phase occours when 
the quark mass is between the two dashed lines {W% + 2We > 0). In addition there is one massive mode inside 
the Aoki phase. 



argued that the Aoki phase and the Sharpe-Singleton scenario potentially both should be realizable 
also in quenched simulations. On the lattice, however, quenched lattice simulations consistently 



observe the Aoki phase Q27| , \2Sj [28| , [3O0, while in unquenched simulations both the Aoki and the 
Sharpe-Singleton scenario are observed [gl], ||, ||, gj, |(J ||, ||, |6|, |^ H H' @> 0I- As 
we now show this puzzle has a natural solution when we consider the problem from the perspective 
of the Wilson Dirac eigenvalues. 



5. Spectrum of the Wilson Dirac operator for W$ = W-j = 

We now derive the spectral density of the Wilson Dirac operator. Let us start with the case 
where = Wj = 0. In this case W% + 2W& > 0, since W% > by the 75-Hermiticity of the Wilson 
Dirac operator, and we have the Aoki phase. 

Because of the y 5 -Hermiticity the eigenvalues of the Wilson Dirac operator are either purely 
real or come in complex conjugate pairs [^2|]. For a scatter plot of the eigenvalues, see for example 
[43]. In order to derive the spectral density of the Wilson Dirac operator in the complex plane 

P cM^*) = (Ls 2 (xy-z)) Nf (5.D 

i 

from Wilson chiral perturbation theory we therefore need to extend the partition function into a 
replicated generating functional with p additional pairs of conjugate quarks with masses z and £*, 

z v Nf+2p (z.XM) = <IW -«)I1& W -zY^TY - *T>- (5-2) 

i i 

The eigenvalue density of the Wilson Dirac operator in the complex plane is then obtained as 

PeJfMt,ih;di) = d r Hmd 2 -logZ^ +2 (£,z*,m;a ; ). (5.3) 
1 p— >o p it- 



6 



Chiral Dynamics With Wilson Fermions 



K. Splittorff 




Re[z}LV 



Figure 4: Left: The microscopic eigenvalue density of the Wilson Dirac operator in the complex plane (in 
the sector with v = for Nf = and a^/W%V = 0.75, W(, = Wj = 0) and Right: the microscopic eigenvalue 
density of the real eigenvalues of the Wilson Dirac operator (in the sectors with v = 1,2,3 for Nf = 
and ay/WgV = 0.2, We = Wj = 0). The exponential fall-off of the eigenvalue density is essential to avoid 
exceptional configurations. 



The trick is now to use the low energy effective theory, ie. Wilson chiral perturbation theory, to 
compute the replicated generating functional and thereby obtain the density. For an introduction to 



the use of the replica method in chiral perturbation theory, see [ }44[ ] . 

The low energy limit of the replicated generating functional is given by Wilson chiral pertur- 
bation theory with the mass matrix 

J{ = diag (m Nf ,z p ,t p ), (5.4) 

where the subscript on the quark masses refers to the number of times the specific mass appears. In 
the mean field limit the dependence on the number, p, of replica flavors is trivial for W$ = Wj = 0, 
and one simply obtains (we use the notation z = x + if) 

P C %= 2 (*A<3 8 ) = 0(8fli - (5.5) 

Note that the eigenvalue density is independent of the imaginary part of the eigenvalue, hence it 
forms a strip along the imaginary axis of width 8<5g. The Aoki -phase sets in when the quark mass 
enters the eigenvalue strip, as in figure |[ 

The microscopic limit of both the eigenvalue density of Dw in the complex plane and on the 
real axis was worked out in [^,[8, 11, 14], see figure Q The exact way in which the real eigenvalues 



are distributed shows at which value of the quark mass the simulations are safe from exceptional 
configurations. 



6. The realization of the Sharpe- Singleton scenario: The solution to the puzzle 
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Figure 5: The realization of the Sharpe-Singleton scenario (Wg + 2W& < 0). The red cross marks the size 
of the quark mass top: the mass is positive and below: the mass is negative. As the quark mass changes 
sign the strip of eigenvalues jumps across the origin and consequently the chiral condensate has a first order 
discontinuity at m = 0, cf. the middle panels. Since the quark mass never reaches the eigenvalue density the 
pions remain massive even at m = 0. See the rightmost figures. Clearly this scenario is only possible in the 
unquenched theory where the eigenvalue density can depend on the quark mass. 



Let us now compute the effect of W(, on the spectral density of Dw in the mean field limit. To 
do so it is useful first to linearize the square in the W(, term of the replicated generating functional 
with a Gaussian integral over a fluctuating mass yg. This leads to [14] 

1 



p C}Nf {z,z ,m;a 6 ,a & ) 



Z^ f (m;d 6 ,d 8 



J [dyt] Z N f (m-y 6 -d & )p^ Nf (z-y6,Z* -y6,m-y 6 -d s ) 



(6.1) 



where we used the notation [dye] = dye/{4\^7l\d(,\)exp(—y^/(l6\d^\)). This expresses the eigen- 
value density of Dw in the complex plane at We / in terms of the eigenvalue density with We = 0. 
At the mean field level for Nf = 2, using Eq. (5.5), we get 

1 



.MF 



where 



ZY F (m;d 6l d & ) 



J ^ 6 ^/ 16 l^lz 2 MF (m-j 6 ;« 8 )e(8^-|i-j 6 |), (6-2) 



2 {m;a 6 ,a&) 



e 2m+\6\a 2 6 \-Aa\ . g -2m+16|flg|-4a| 



+e(8(ai + 2a£)-Hy 



(6.3) 



m 2 /8(al-2\dH)+4aj 
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This saddle point structure of the two flavor partition function leads to 

1 



„ MF / *■ \ 



Z^ F (m;a 6 ,a 8 ; 



(6.4) 



X j e 2 ™+ l 6\al\-4al 



d(Sal-\x+\6\a 6 \ 2 
+6(S(al+2dl)-\m\)6 [Sdj 



x — 1 6 1 «6 1 2 1 



x + 



2\d(,\ m 



(a| — 2|a 6 | 2 ) 



> 2 /8(«1-2|«1l)+4a| 



Since we seek to understand the realization of the Sharpe-Singleton scenario let us consider Wg + 
2W(, < 0. The density in this case is shown in figure |5[ As the quark mass changes sign the strip of 
eigenvalues jump from one side of the origin to the opposite side and this causes the 1st order jump 
of the chiral condensate. This explains how the Sharpe-Singleton scenario is realized in terms of 
the strong mass dependence of the unquenched eigenvalue density. 

In the quenched case on the contrary the eigenvalue density is independent of the quark mass 
since the fermion determinant is absent in the measure. The quark mass in the quenched case, 
therefore, necessarily passes through the eigenvalue density and this leads to the standard Aoki 
phase, independent of the sign ofWg + 2W(>. This explains why the Aoki phase is observed con- 
sistently in quenched simulations and thus solves the puzzle of section f| One can also compute 
the quenched chiral condensate directly from the graded version of the partition function in Wilson 
chiral perturbation theory and reach the same conclusion. This is illustrated in figure H. The anal- 



ysis of [26] reached the conclusion that the Sharpe-Sharpe-Singleton scenario could be realized 
in quenched simulations because the bound on the signs of the Wi's due to 75-Hermiticity had not 
been understood at the time. 

The meta-stable region surrounding the 1st order Sharpe-Singleton phase transition (see figure 
|2|) can also be observed in the eigenvalue density. For < mL < \ Wg +2W(,\a 2 and W& + 2Wg < the 
meta-stable minimum of Z^ F introduces a local minimum of the % integral in Eq. (6.4). This local 



minimum offers the possibility for the strip of eigenvalues to temporarily jump to the opposite side 
of the origin thus inducing a fluctuating sign of the chiral condensate. Such fluctuations become 
more frequent as the meta-stable and global minimum of Z^ F become almost degenerate, ie. as m 
approaches zero. Thus the meta-stability of the eigenvalue density is completely consistent with 



the meta-stability of the chiral condensate obtained directly from Eq. ( p.4| ). 

It is also instructive to consider the distance from the quark mass to the edge of the strip of 



eigenvalues of Dw- From Eq. (6.4) we see that this distance is given by 



\m\ -8(Wg + 2W 6 )a 2 /L. 



(6.5) 



This is exactly the combination which enters the right hand side of the order-a 2 corrected GOR 



relation, cf. Eq. ( |3.5| ). The distance from the quark mass to the edge of the strip of eigenvalues of 
Dw can therefore be thought of as the effective quark mass which enters the standard continuum 
form of the GOR relation. 

In the Sharpe-Singleton scenario, ie. for Wg + 2W$ < 0, the quark mass never reaches the strip 
of eigenvalues of Dy/. The minimal distance between the quark mass and the eigenvalues is given 
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Figure 6: The Sharpe-Singleton first order phase transition is only realized in the unquenched theory. 
Here we show the mass dependence of the microscopic chiral condensate for ag = 1 and &(, = 0,0.5/ and ; 
corresponding to W% + 2W 6 > 0, Wg + 2W 6 = and Wg + 2W 6 < 0, respectively. Left Nf = 2: The 1st order 
Sharpe-Singleton phase transition develops with the increasingly negative value of We- Right Nf = 0: The 
smooth behavior of the chiral condensate as a function of the quark mass characteristic of the Aoki phase is 
present in the quenched theory even when W% + 2W(, < 0. 

by —8(^8 + 2W^)a 2 and hence the smallest value which the pion mass can reach is 

m n = ^-16(W 8 + 2W 6 )a 2 /Fl (6.6) 

This is of course exactly the same minimal value of the pion mass found in [g]. We now understand 
how this minimum is linked to the behavior of the eigenvalue density in the unquenched theory. 



7. Spectrum of the Hermitian Wilson Dirac operator 

The 75-Hermiticity of the Wilson Dirac operator makes it natural to introduce the Hermitian 
Wilson Dirac operator, D5, defined by 

D 5 = y 5 (D w + m) . (7.1) 



There is a close analogue of the Banks Casher relation p5[ ] for the spectrum of D5: The lattice 
artifacts induce eigenvalues of D5 with magnitude less than \m\ and when these build up a density 



at zero the Aoki phase is reached Q46|]. It is therefore essential to understand analytically the 
dependence of the smallest eigenvalues of the Hermitian Wilson Dirac operator in order to study 
chiral dynamics in simulations with Wilson fermions. 

The first to consider the spectral density of D5 from the perspective of Wilson chiral perturba- 
tion theory was Sharpe 0. By now a detailed analytic understanding of the effects the discretiza- 
tion eiTors have on the spectrum of D5 have been obtained [||, ||, ^, 10, 11, O, O, 14, [l5|, [16|, 0|. 
A plot of the spectral density is shown in figure [7|. The plot shows the density in the sector with 
V = 1 and the index peak at m is clearly visible. Also the precise way the eigenvalues intrude into 
the region between — m and m can be observed. 
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Figure 7: Red histogram: The spectral density of the Hermitian Wilson Dirac operator on quenched 
20 4 lattices in the sector with v = 1 . The eigenvalues have been rescaled with TV /a= 220. See [ |50| ] for 
details about the simulations. Blue curve: The microscopic spectral density obtained from Wilson chiral 
perturbation theory with Nf = 0, V = 1, m = 2.3 and d(, = 0.06. Also shown is the black curve obtained 
from ordinary chiral perturbation theory with a = 0. The prime effect of the lattice artifacts for ha\ <C 1 
is to smear the index-peak. This leads to the a / \fV scaling of the width of the smallest eigenvalue of the 
Hermitian Wilson Dirac operator reported in [|Ti| |i"9| |. 




-20 -10 10 20 



Figure 8: The accumulated single eigenvalue distributions of the lowest eigenvalues of the Hermitian 
Wilson Dirac operator. Plot adopted from [p^l. 



The predictions also offer a new way to measure the low energy constants of Wilson chiral 
perturbation theory [17, 48, 49, 50]: The histogram in figure ^ displays lattice data and the match 
of the analytic prediction to the lattice data fixes the values of the low energy constants. We refer 
to Q for details. 

One can also obtain the low energy constants from a fit of the analytic predictions for the 
spectral density of Dw to lattice data. However, since D5 is Hermitian it is somewhat easier to 



determine this spectrum numerically in lattice simulations. See also pl\ ] for an alternative way to 
measure the W,'s. 



8. Scaling with aj \jV 



In the simulations of I 



19, 20(1 it was found that the width of the distribution of the smallest 



eigenvalue of D5 scales with a/y/V. At first sight this may sound disturbing, since the Banks- 
Casher relations [45, ^6] suggests that the smallest eigenvalues should scale with 1 /V. Hence one 
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Figure 9: The accumulated eigenvalue density the Hermitian Wilson Dirac operator including all sectors. 
Left: quenched lattice simulation from [p2[]. Right: example of the analytic prediction from Wilson chiral 
perturbation theory in the e regime with Nf = 0. The parameters chosen for this example are mLV = 3, 
ag = 0.2 and W 6 = W 7 = 0. 



might fear that the a/\/V -scaling is a sign of uncontrolled lattice artifacts, such as local defects not 
described by Wilson chiral perturbation theory. On the contrary, the a/y/V scaling is a good sign: 
In the limit where d{ <C 1 the smallest eigenvalues of D5 should in fact scale precisely with a / \fV 
[§, §]. For di <C 1 the primary influence of the lattice artifacts is to smear the topological mode 
into a peak around m of width a / \/V , as shown in figure This fact, obtained from Wilson chiral 
perturbation theory, has been demonstrated explicitly in lattice quenched simulations [ff8|, ^9j, [50| ]. 
See figures [7] and |[ Note that these lattice simulations were separated into sectors with fixed 
index. The effect of the discretization error on the spectrum of D5 can also be studied without 
the separation into sectors. In figure the accumulated eigenvalue density including all sectors is 



shown. The left hand plot is lattice data obtained in [ |52| ] while the right hand plot gives an example 
of the analytic results. It would be most interesting to preform a systematic fit to the data including 



the leading order corrections to the slope at large A obtained in [17]. 



9. Alternative vacuum in the Aoki phase for N/ = 2 



In [53] it was proposed that the additional condensate (iujsu + idjsd) will take a nonzero 
expectation value inside the Aoki phase along with the standard Aoki condensate (iujsu — idysd). 
This suggestion is particularly interesting since the presence of the additional condensate is in direct 
contradiction with the results of Wilson chiral perturbation theory [|54]]. The additional condensate 



(iuy^u + idy^d) is zero in the £-regime of Wilson chiral perturbation theory, Eq. ( |2.3[ ), because 
it is the v.e.v. of Tr(£/ — £/') which vanishes in SU(2) where TrU = TrU^. With fixed index, 
however, the group average is extended to U(2) and Tr(U — W) does not vanish trivially. The 
behavior of {iuy^u + idy=,d) at fixed index can be examined by explicit evaluation of U (2) integral 
Eq. ( |2.5[ ) and differentiation w.r.t. £. One can also evaluate the squared condensate ((iujsu + 
idysd) 2 ) at zero external sources as suggested in [55]. Again this is zero in SU(2) and nonzero 



at fixed index. As should of course be true, one recovers the vanishing value of the additional 
condensates within Wilson chiral perturbation theory, after a careful summation over all sectors 
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with the appropriate partition functions as weights. In simulations, however, it can be very tricky 
to control this sum over the sectors. Even a small error in the sampling of the sectors can lead 
to a substantial error in the value of the additional condensates after summation over the sectors. 
The summation over the sectors become particularly delicate in the Aoki region if the low energy 
constant W-j is non-vanishing. The reason is that the term proportional to W-j is tantamount to a 
Gaussian fluctuating source for iuy^u + idy^d . In other words, Wilson chiral perturbation theory 
predicts that it is essential to sample all the sectors with the correct weight in order to determine if 
the additional condensate is present. 



10. Conclusion 

The approach to chiral dynamics with Wilsons fermions from the perspective of the eigen- 
values of the Wilson Dirac operator has lead to several new insights. On the macroscopic level 
it becomes clear that the Sharpe-Singleton scenario can only be realized in unquenched theories 
and on the microscopic scale we learn that the a / y/V scaling of the smallest eigenvalues of the 
Hermitian Wilson Dirac operator is a sign that the simulation is very close to the continuum. On a 
more technical but equally important level we now understand that the signs of the additional low 
energy constants in Wilson chiral perturbation theory are fixed. Moreover, the new analytic results 
offers a practical way to measure the additional low energy constants in lattice simulations. These 
combined new insights can be most useful when we seek to minimize discretization errors in lattice 
simulations with Wilson fermions. For example, to minimize + it is useful to know that the 
two contributions have opposite signs. 

A control of the lattice artifacts becomes particularly important in studies of QCD with a larger 
number of flavors: In a mean field approach the coefficient in front of W(, and Wj scales with Nj 
compared to the linear scaling of the W% term. This suggest that simulations with a larger number 
of flavors should be more likely to end up in the Sharpe-Singleton scenario. Hence one is less likely 
to be able to exploit the massless modes at the boundary of the Aoki phase to mimic the chiral limit 
for larger number of flavors. 

Many aspects discussed in this review of chiral dynamics with Wilson fermions has a direct 
analogue for staggered fermions. It would be interesting to pursue this further. The necessary tools 
are already available and first results have appeared [p7|]. The spectral properties of the Dirac 



operator could also cast additional light on simulation with mixed actions, see eg. [58]. 
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